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Abstract 

The most general spin structures of the spin-1/2 baryon and spinless meson production operator 
for both photon and nucleon induced reactions are derived from the partial-wave expansions of 
these reaction amplitudes. The present method provides the coefficients multiplying each spin 
operator in terms of the partial-wave matrix elements. The result should be useful in studies of 
these reactions based on partial-wave analyses, especially, when spin observables are considered. 
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In the present work, we derive the most general spin structure of the reaction amplitude 
for both positive and negative parity spin- 1/2 baryon production in nucleon-nucleon (NN) 
collisions and also in photon induced processes on nucleons. Knowing the spin structure 
of the transition operator is of particular importance in analyses of spin observables. The 
method used here to extract the spin structure is a generalization of the partial-wave ex- 
pansion of the NN amplitude following Ref. Q]. It is quite general and, in principle can be 
applied to any reaction process in a systematic way Usually, the structure of a transition op- 
erator is derived based solely on symmetry principles. The usefulness of the present method 
is that it also provides explicit formulae for the coefficients multiplying each spin struc- 
ture in terms of the partial-wave matrix elements and this should be particularly useful in 
model-independent analyses based on the partial-wave expansion of the reaction amplitude. 

As an example of application of the present method, we investigate the possibility of de- 
termining the parity of a narrow resonance in the photon as well as in the nucleon induced 
reactions. The search for new resonances, especially the so-called missing resonances and 
exotic resonances, are receiving increased attention |2|. Apart from establishing their exis- 
tence, the determination of their basic properties are of extreme importance. Among these 
properties, the parity is of particular interest in connection with the substructure of these 
resonances ^J. However, it is often the case that no theoretical predictions can provide a 
conclusive result for the parity and other basic properties. A recent example of this situation 
is provided by the pentaquark @ + . The existence of this exotic baryon and the determi- 
nation of its spin and parity quantum numbers have been under an intensive investigation 

n 

both experimentally and theoretically over the past couple of years j^J . It is, therefore, 
very important to find a way of determining these properties in a model independent way. 

In the past few years a considerable amount of data for meson production in NN collisions 
have been obtained (see Ref. p for a review). In particular there are data, not only for cross 
sections but, in the case of pion production in NN — > NNtt, a large set of spin observable 
data [7]. In addition, the data bases for the production of other mesons, such as r\ and 7/, 



are growing rapidly 



. Also, reactions involving the production of particles containing 



strange quarks such as the NN — > YNK reaction where Y stands for a hyperon, are receiving 



1 The existence of the G + is still controversial. While this resonance has been seen in a number of experi- 
ments, it has not been observed in a similar number of experiments, including a very recent experiment 
with much better statistics. 
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increased attention [y, |9|. With these (high-precision) data becoming available, there is a 
demand for more thorough and detailed theoretical analyses of these reactions. Therefore, in 
the present work, we also derive the most general spin-isospin structure of the NN — > NB'M 
transition operator, where M stands for a spinless (scalar or pseudoscalar) meson and, B', 
for a spin- 1/2 baryon with positive-parity. The spin structure of the NN — > NNir transition 
operator has been derived in Ref.|ly] and applied to neutral pion production. Here, we also 
consider charged meson production. 

The present paper is organized as follows: In section I, the general spin structure of 
the transition operator is derived for the reaction 7 + iV — > M + B, where B stands for a 
spin-1/2 and either a positive or negative parity baryon. In section II the results derived in 
section I are illustrated by applying them to the near-threshold kinematic regime. Sections 
III and IV are devoted to the derivation of the most general spin structure of the reaction 
iV + iV — > B' + B, where B' stands for a spin-1/2 and positive-parity baryon, and to its 
application near-threshold, respectively. The reaction N + N^M + B' + N is considered 
in sections V and VI. A summary is given in section VII. Appendices A-D contain some 
details of the derivation of the spin structure of the transition operators. 



I. THE REACTION 7 + N M + B. 

We start by making a partial-wave expansion of the 7 + N — > M + B reaction amplitude. 
Here M stands for a pseudoscalar meson and B, a spin-1/2 baryon. We, then, have 

<^m'\M(q,k)\SM s >= £ i L - L ' (SM S LM L \ JMj){±m'L'M L ,\ JMj)M J L ? L (q, k) 

x Y L , MLl (q)Y* ML (k) , (1) 

where S, L, J stand for the total spin, total orbital angular momentum and the total angular 
momentum, respectively, of the initial 7A state. Ms, Ml and Mj denote the corresponding 
projection quantum numbers. The primed quantities stand for the corresponding quantum 
numbers of the final MB state. The summation runs over all quantum numbers not specified 
in the left-hand side (l.h.s.) of Eq.(£Q). k and q denote the relative momenta of the two 
particles in the initial and final states, respectively. The partial-wave expansion given above 
is, of course, related to the more commonly used electric and magnetic multipole expansion. 
For the present analysis, however, it is convenient to use the above expansion. 
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Eq.flf can be inverted to solve for the partial-wave matrix element Ml, J L (q, k). We have 
Mli(q,k) = Y,i L '~ L (l m ' L ' M L>\JMj)(SM s LO\JMj) 



^7T 



2 



2 J + 1 V 4vr 



2L + 1 r +1 „ ,„ u „, ,„ „ x i 



x T7^\/^— / d(cos(0))3£ Mr ,(0,O) < -m'|M(g,fc)|SM s > , (2) 



-i 



where, without loss of generality, k is chosen along the z-axis and q in the xz-plane; cos(#) = 
q ■ k. The summation is over all quantum numbers not specified in the l.h.s. of the equation. 
The most general spin structure of the transition operator can be extracted from Eq.(£Q) 

as 

M(q,k) = \\ m ' >< \rn'\M{q,k)\SM s >< SM S \ . (3) 

SM s m' 

Inserting Eq.((TJ) into Eq.® and re-coupling gives 



M(?,^=X; <W (-)" J " i W a ^i(9.*)E|^ i [)\B s ®A h r.[Y L {k)®Y L ,{q)}\ (4) 



2 a 



where we have used the notations Bsm s = ( _ ) < S — M s \, Ai m , = \\m' >, and 



[J] = \/2J + 1. The outer summation is over the quantum numbers S, L, L', and J. In the 
above equation S is either 1/2 or 3/2 so that a takes the values 0, 1, and 2, and denotes 
the rank of the corresponding tensor. In the above equation it should be understood that 
the matrix elements of the meson creation and photon annihilation operators have already 
been taken. 

We now expand \Bs <E> A for each tensor of rank a, in terms of the complete set of 
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available spin operators in the problem, i.e., the photon polarization vector e and the Pauli 
spin matrix a together with the identity matrix. The result is 

{B s ®Aif = 4j*-iT, 

[B s ®A,f = ^®ef, (5) 

where the numerical factors are uniquely determined such that the spin matrix elements of 
the right-hand side (r.h.s.) in the above equations equal the corresponding matrix elements 
of the l.h.s. 



4 



What we have done so far applies to either a negative or positive parity baryon B. Total 
parity conservation demands that (— ) L+L ' = +1 and (— ) L+L ' = —1 in the case of a positive 
and negative parity B, respectively. This leads to distinct spin structures of the transition 
operator for positive and negative parity as we shall show below. Hereafter, the superscript 
± on any quantity stands for the positive (+) or negative (— ) parity of B. 

A. Positive parity case 

For a positive parity baryon B, choosing the quantization axis z along k, the quantity 
[Yz,(/c) ® YL'{q)] a can be expressed without loss of generality as 

[Y L (k)®Y L ,(q)}° = H l[ ^Pl (q ■ k)6 L>L , , 

[Y L {k)®Y L ,{q)Y = l Hf 

[Y L (k) ® Y L ,(q)} 2 = a VL [q <g> q} 2 + b UL [k ® k} 2 + c VL [k <g> q} 2 . (6) 

The structure in the above equation is dictated by total parity conservation. Pl(PI) is the 
ordinary (associated) Legendre function. n 2 = (k x q)/\k x q\. The coefficients ajyi, b^i 
and cl'l are derived explicitly in Appendix A. 
Inserting Eqs. (j5l6p into Eq.(0J) we have 

M + (g, k) = T\B ■ t + iT^e ■ n 2 + Tzip x e) • h 2 

+ T±\o ® e\ 2 • [q <g> g] 2 + JS[^ ® ef • [fc <g> A;] 2 + T^\a ® t\ 2 ■ [k <g> q} 2 , (7) 
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where 

^ 3 = 7!(wtt)^' 

The summations are over all quantum numbers appearing explicitly in the r.h.s. of the 
equalities. 

The above result is the most general form of the spin structure of the transition operator 
consistent with symmetry principles. The coefficients ft are functions of the energy of the 
system and scattering angle 9 of the meson M relative to the photon direction. We note 
that the transition operator in Eq.((7|) has six independent spin operators and is valid for 
both real and virtual photons. Using the identity 

3 

3[ a® e ] 2 ■ [ a <g> b f = -[ a ■ ae ■ b + a ■ be ■ a ] - (a ■ b)a ■ e , (9) 
where a and b stand for arbitrary unit vectors, Eq.tJTJ) can be rewritten as 

M + (q, k) = Fi<? • e + iF 2 e ■ h 2 + F 3 a -ke-q + F 4 a -qe-q + F 5 a -ke-k + F 6 a -qe-k , (10) 
where 

1 r .1 
ft + ft + (g . k)ft 



ft — ft , F 3 — -j -ft 3 + -ft 

\k x q\ * 



F A = ft, F 5 = ft, F 6 = ~-J—ft + \ft ■ (11) 

\k x q\ * 

It should be noted that Eq. ()10|) is equivalent to that of Ref . [ill] . In fact, apart from an 
irrelevant overall factor of i, the coefficients fj = jF y< - ±,0 \ (j = 1, 6) in Ref.|lj| are related 



to those in Eq.ffTU]) by 

Fx = fx - {q ■ k)h , F 2 = \kx q\f 2 , F 3 = (/ 2 + / 3 ) , 

F 4 = / 4 , F5 = -A-(ff.fc)/3-^/ B , F 6 = -(q ■ k)U - £/ 6 . (12) 

fc Ko 

In the case of a real photon (photoproduction amplitude), the number of independent 
spin operators in Eq^] reduces to four due to the transversality condition, i. e., the terms 
_F 5 and _F 6 are absent. EgjlUI then becomes equivalent to that of Ref.Q]. The coefficients 
fj,(J = 1,...,4) in Ref.(l2j are related to the corresponding coefficients Fj in Eq. ffTUj) as 
given by Eq. (fT2| with k 2 = e ■ k = 0. 
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12j is that in the present work 



One difference between our work [Eq. (Jl(J|) ] and Refs. 
these coefficients are related explicitly [Eq.(|HJ)] to the partial- wave matrix elements intro- 
duced in Eq.([T]). 

B. Negative parity case 

For a negative parity baryon B, choosing the quantization axis z along k, [Yl(A;) ®Yl> (q)] a 
can be expressed as 

[Y L (k)®Y L ,{q)f = 0, 
[Y L {k)®Y L /^ :LL ' 



i 7T 

\Y L {k) Y v {q)} 2 = a! VL [k ® h 2 } 2 + b' VL [q ® h 2 } 2 , (13) 

where hi = \{k x q) x k]/\k x q\. The coefficients a' L , L and b' L , L are calculated explicitly in 
Appendix A. 

Inserting Eqs. (j5ll3J) into Eq.(JU) we have 

M~(q, k) = iQie- hi + Q 2 a ■ (e x hi) + Q 3 a ■ (e x k) 

+ Qi\a ®e\ 2 -[k® h 2 f + G 5 [a <g> e] 2 • [q h 2 } 2 + %Q % e- k , (14) 



LI , u + 1 A LUL'l\ll)P l L ,(k ■ q)hi + (L0L'0\lQ)P L ,(k • q)k 
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where 



<h = ^ ( ^1 v^(-r^r , [LL ' ] crop \ s L J 

x M J L ? L (q,k)Pb(k-q) , 




£3 = ^= & 



/6 



L' i 2 



I' A 2 

S 6 = i- (i±^) j:^(-)-*[LLl[J]-(LQL'0|10) {* J J 

X M J L ? L (q,k)P L ,Ck-q) . (15) 

The summations are over all quantum numbers appearing explicitly in the r.h.s. of the 
equalities. 

The above result is the most general form of the spin structure of the transition operator 
for a negative parity baryon B. As in the positive parity case, we note that the amplitude 
in Eq. has six independent spin operators and is valid for both real and virtual photons. 
Using Eq.Q, it can be rewritten as 

M~(q, k) = id e ■ q + G 2 a ■ (e x q) + G 3 a ■ (ex k) 

+ Gi [a ■ ke ■ hi + a ■ h 2 e- k) + G§ [a ■ qe- h 2 + & ■ h 2 e- q] + iG% <?• k , (16) 

where 

Gi = 7r^& > G 2 = 7^—-G2 , G 3 = G 3 -(k- q)G 2 , 
\k x q\ \k x q\ 

G A = l -g A , G 5 = ~g 5 , G, = g Q -{k-q)G 2 . (17) 

Quite recently, Zhao and Al-Khalili k| have also given the spin structure of the pho- 
toproduction amplitude in connection to the reaction — > KQ + for the case of negative 

8 



parity 6 + . The structure given in Eq. flTSj) with e • k = is equivalent to that of Eq.(18) in 



Ref. 13[, except for the term a-h 2 e-q which has not been included in Ref. 



13j | on the grounds 



that it is a higher-order contribution. However, this term and the a ■ qe- h 2 term contribute 
with the same coefficient We also note that the recent model-independent analysis of 
the + photoproduction [lj] has been based on the present results and, in particular, on 
Eqs.flHIIini) with e • k = . 



II. APPLICATION: NEAR THRESHOLD AMPLITUDE IN jN MB. 

As an application of the present results, we consider the reaction jN — > MB in the near- 
thresholds kinematics. The complete transition amplitude is given by Eq. ljlOj) for a positive 
parity baryon B and by Eq. ljlfij) for a negative baryon B, respectively, with e ■ k = 0. In 
the near-threshold energy region, the final MB is mainly in relative S and P waves. Then, 
considering only L' = 0, 1, there are seven partial- wave amplitudes. For positive parity B, 
they are 

1 S 1 - S x , 3 D t - S, , ^Pt - P 1 , Vp 3 ^ P 3 , 3 F 3 ^ P 3 , (18) 

where we have used the notation 2S L 2 j — > L' 2J . For these amplitudes, the coefficients 
Fi, (i = 1, 4) in Eq. ffTUj) exhibit the following angular and energy (due to q L ') dependences 



F 1 = A + (|) A x cos(^) , F 2 = (|) ^ sin(0) , F 3 = (|) d , F 4 = . (19) 

In the above equation, Aj/, -B^/ and Cjy denote the linear combinations of the partial- wave 
matrix elements Mjf L resulting from Eqs. f|8lll() for those states specified in Eq.(|18|) having 
orbital angular momentum L' . They are given explicitly in Appendix B. Note that the factor 
(q/A) L ' contained in Mjf L due to the centrifugal barrier has been displayed explicitly in the 
above equation. A is a typical scale of the problem which may be taken to be the four- 
momentum transfer at threshold, A 2 ~ — t = [tub — m 2 N /{rriB + tom)]™m, where tub, thm 
and m^v denote the masses of the baryon B, meson M and nucleon, respectively. Therefore, 
near threshold, higher partial- wave contributions will be suppressed by the factor {q/K) L ' 
if heavy particles are produced in the final state. Moreover, for short-range processes, the 
coefficients Ay, B L i and Cy are nearly constant independent of energy. 
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Analogously, for a negative parity baryon P, the possible partial-wave amplitudes are 



3,1 



> L P X -> S x , l Si -> Pi , 3 L>x -> Pi , d S 3 -> P 3 , ^Dg -> P 



3 ■ 



(20) 



With these amplitudes, the coefficients G% in Eq. ()16|) exhibit the following angular and 
energy (due to q L ') dependences 



G 3 



A 



4 



A 



p; cos(^) 



G 2 

G4 : 



V2- 1 

(|)Gisin(0), G 5 = 0. (21) 



In the above equation, A' L ,, B' L , and C^/ denote the linear combinations of the partial- wave 
matrix elements Mjf L resulting from Eqs. (jl5ll7j) for those states specified in Eq. (|20|) having 
orbital angular momentum V . They are given explicitly in Appendix B. For a short range 
process, A' L ,, B' L , and C' L , are nearly constant independent of energy 

Following Ref. Q| we now introduce e± = y and eji = x denoting the photon polarization 
perpendicular to and lying in the reaction plane (xz-plane), respectively. Recall that the 
reaction plane is defined as the plane containing the vectors k (in the +z-direction) and q 
and that k x q is along the +?/-direction, in which case, fi\ = x and fi2 = y. Then, from 

Eas. ffTnnn 



M + 



A + ( ± ) A x cos{6) 



<y,, 



l) Bx M&) , 



A +[j)Ax cos(0) 



cr x + 



Similarly, from Eqs. (jl6l21|) 



M-ll 



/ (fj ^cos(0) 
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JJCisin(0) a 2 



|) CI sin(£)a 2 

q 



B'o+[j) B[ cos(#)j a, + * j ^ sin(0) . 



(22) 



(23) 



In the above equation, B[ = B[ + A[ and C[ = C[ - A[, with A[ = ^ (75^7) ^'i- 

Following Ref.Q], any observable in the 7 A — > MB reaction can be readily calculated 
from Eqs. ()22I23J) for both positive and negative parity baryons P. It is then straightforward 



to show that no observable in this reaction is able to distinguish between positive and neg- 
ative parity P, unless one measures the polarization of P. In particular, neither the energy 
dependence nor the angular distribution exhibits features sufficiently distinct to determine 
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the parity of the baryon B unambiguously. As has been shown in Ref. Q| in connection 
with NN — > YQ + , the situation is quite different in the NN — > B'B reaction. We shall 
discuss this latter reaction in the following two sections. 



III. THE REACTION N + N -» B' + B. 

We now consider the process N + N —>■ B' + B, where B' stands for a positive parity spin- 
1/2 baryon and B, either a positive or negative parity spin-1/2 baryon. The partial- wave 
expansion of the corresponding reaction amplitude is 

<S'M s ,\M(f',0)\SM s >= £ i L - L '{S'M s ,L'M v \JMj){SM s LM L \JMj)Mi;, s L J (p',p) 

x Y l , Ml ,(p')YZ Ml (p) , (24) 

where S, L, J stand for the total spin, total orbital angular momentum and the total angular 
momentum, respectively, of the initial NN state. Ms, Ml and Mj denote the corresponding 
projection quantum numbers. The primed quantities stand for the corresponding quantum 
numbers of the final B'B state. The summation runs over all quantum numbers not specified 
in the l.h.s. of Eq. (j24)l . p and p ' denote the relative momenta of the two particles in the 
initial and final states, respectively. We note that, in Eq. (j24j) . apart from the restrictions on 
the quantum numbers encoded in the geometrical factors, total parity conservation imposes 
that (— ) L+L ' = +1 and (— ) L+L ' = —1 for positive and negative parity B, respectively. 

Eq. (j24J) can be inverted to solve for the partial-wave matrix element M£,'£ J (p',p). We 
have 

Ml'l J (p',p) = ^i l '- l {S'M s ,L'M l ,\JMj)(SM s LQ\JMj 



57T 



2J+ 1 



or -i_ i r +i 

x d(cos(9))Yl /ML X9,0)<S'M s/ \M(p',p)\SMs> , (25) 

where, without loss of generality, the z-axis is chosen along p and p ' in the res-plane; 
cos(6 l ) — pi -p. The summation is over all quantum numbers not specified in the l.h.s. of the 
equation. 

The most general spin structure of the transition operator amplitude can be extracted 
from Eq. (J2U) as 

M(p',p) = \ S ' M S' >< S'M s >\M(p',p)\SMs >< SM S \ . (26) 

S'SM S M S , 
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Inserting Eq . ()24p into Eq. ()26|) and re-coupling gives 

M(P',P) = "£i L - L 'n J+s '[j] 2 Mi; s L j ( P ',p)j2 



S L J 
V S' a 

x [A g , ® B s f ■ [Y L ,(p') ® Y L (p)] a , (27) 



where we have used the notations B SMs = (— ) s ~ M s < S — M$\ and As>m s , = \S'M S / >. 

We now expand [Agi ®> B$] a , for each tensor of rank a, in terms of the complete set of 
available spin operators in the problem, i.e., the Pauli spin matrices d\ and <? 2 , corresponding 
to the interacting particles 1 and 2, together with the identity matrix. Then, a takes the 
values 0, 1, and 2, and denotes the rank of the corresponding (spin) tensor. There are six 
cases to be considered: 



S = S' = 0, a = : [A s > <g> B s }° = |00 >< 00 1 = P< 



5=0 



S = S ' = l,a = 0: [A S ,®B S ]° = -^=J2\ 1M s ><1M s\ = -/ 



S = S ' = l,a = l: [A s > <S> Bs] 1 = ^(^ + ^2) - 

5 = 0, 5' = 1, a = 1 : [A S / ® Ps] 1 = J (^1 - a 2 ) Ps=o 



5 = 1, 5" = 0, a = 1 : (8) Ps] 1 = -1 (^ - <r 2 ) Ps=i , 

S- = 5" = 1, a = 2 : W 5 / <g> P 5 ] 2 = -Lfo ® ^ 2 ] 2 , (28) 

where Ps stands for the spin projection operator onto the (initial) spin singlet and triplet 
states as S = and 1, respectively. In terms of the Pauli spin matrices we have Ps=o — 
(l-a 1 -a 2 )/4andP 5=1 = (3 + ^ -a 2 )/4. Also, {a 1 -a 2 )P s /2 = [(a x -a 2 ) + {-) s i{ffi x a 2 )]/4. 

In the following we shall consider the case of a positive and a negative parity baryon B 
separately. 



A. Positive parity case 

For a positive parity baryon B, the quantity [Yl>(p') ® Yk(j3)] a in Eq. (j27|) can be read off 
from Eq.©. Note that \Y L >{p') ® Y L (p)] a = (-) a [Y L (p) <g> Yl'(p')] q for the positive parity 
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case. Inserting this and Eq. (|2*5j) into Eq. (|27jl . we have 

M + (p\p) = V 1 P S=0 + V 2 P S=1 

+ iVa(B x + ^2) • fi 2 + iD^Bx - a 2 ) ■ n 2 P s=0 + %V^{a x - &*) ■ n 2 P s=1 

+ P 6 [ai <g> a 2 ] 2 ■ [p <g> p] 2 + ©7^! <g> a 2 ] 2 • [p' g) p'] 2 + Pg^i ® ^a] 2 • [p ® pT 



(29) 



where 



47T 3 z — ' 



X> fi 



©7 



|E^ X 'H J+1 [^] 2 M^(p',p)a^ 
2 ^ L 1 2 



V 1 2 



1 ^ L - L \-) J+1 \Jf M}}l(p>,p)c L , L , (30) 

2 ^ L' 12' 



The coefficients az/L, and cx/z, are given by Eqs. f|77|78^ with the replacements k — > p 
and q —> p'; these same replacements are needed to calculate h 2 in Eq.(j2HJ). 

Eq.(|29|) is the most general spin structure of the transition operator for a positive parity 
baryon B, consistent with symmetry principles. It contains eight independent spin struc- 
tures. The first two terms are the central spin singlet and triplet interactions, respectively. 
The third term is the spin-orbit interaction. The fourth and fifth terms describe the spin 
singlet— ^triplet and triplet— ^singlet transitions, respectively. The last three terms are the 
tensor interactions of rank 2. Apart from the fourth and fifth terms, all the other terms 
conserve the total spin in the transition. 
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It should be mentioned that in the case of identical particles, i.e., NN — > NN, the 
structure (<7i — o-i) is not allowed in the scattering amplitude, which reduces the total number 
of independent spin structures in Eq. (j2T?j) to six [l]. Furthermore, for elastic scattering, 
T>q = T> 7 , as a consequence of time reversal invar iance 

For the purpose of calculating the observables directly from the transition operator of 
Eq.(|29|). it is convenient to re-express it in the form 

9 3 

M + = Y, E KVn(lK<(2) , (31) 

A=l n,n'=Q 

where <Jo{i) = 1, cri(i) = cr x {i), etc., for i-th nucleon. The coefficients Ml are linear 
combinations of the coefficients appearing in Eq. (j30|) . Explicitly, we have 

1 



M l nn , = - [3V 2 + V,] 6 nfi 6. 



n'O 



M 3 nn , 



V 3 + -{V 4 + V 5 ) 



^2n(l — 8n,0)8n',0 , 



ML, 
ML, 



ML, 



ML' 



— [X> 5 — X> 4 ] £ nri 'fcn2fc(l — ^n,o)(l — ^n',o) 

T>sp n p n >(l - ^n,o)(l ~ #n',o) , 
^7P^(l-5 n ,o)(l-5 <0 ) , 
-T>8pnP' n ,(l ~ 5n,o)(l - <*n',o) , 
^V 8 p n p n ,(l - <y n , )(l - £n',o) , 



(32) 



where we have used the notation a n for the n-th component of an arbitrary unit vector a; 
e nn 'fc denotes the antisymmetric Levi-Civita tensor and A is an index for the type of term 
(operator) being considered. In obtaining Eq. (j32j) . we have made use of the identity 



3[ (?i ® &2 ] 2 ■ [ o ® S ] 2 = -[ ai • ac? 2 • S + CTi • 6a 2 • a ] — (a • 6)<Ti • e?2 



(33) 



where a and 6 denote arbitrary unit vectors. 

With the transition operator in the form of Eq. (J31|) . it is straightforward to calculate any 
observable of interest. Of course, it can also be expressed in terms of the matrix elements 
of total spin of Eq.(|24|): for convenience some observables are given in Appendix D in both 
representations. 
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B. Negative parity case 



For a negative parity baryon B, the quantity [Y/yfj?') (8>lL(p)] a in Eq. (}2Tj) can be read off 
from Eq.(JT3J). For this case [Y L r(fi) ® Y L (p)) a = {-) a+1 \Y L (p) <g> Y v {p)] a . Inserting this and 
Eq.flHI) into En.(f27jl. we have 

M~{p',p) = i{H 1 n 1 + H 2 p)-(a l + a 2 ) 

+ «(?4ni + 7^4j5) ■ (^l - o 2 )Ps=o + i(H 5 ni + 7i 6 p) ■ {B x - a 2 )P s =i 

+ H 7 [ai ® a 2 ] 2 ■ [p ® n 2 ] 2 + ft 8 [^i ® ^] 2 • [p ® n 2 ] 2 , (34) 

where 

^ = ^^E^ L ' +1 (-) J [ LL, n j ] 2 ( L0L,0 i 10 )|^ j ^|m^(p',p)p^(p'-p), 

Ws = ^ E* L ~ £ ' +1 (-) £/+1 -Jp= 1} (^oriiii)Mfff 0/,^^ • p) , 

2 ^ I^L' 1 2 J 

^ = iE^'HW I 1 ^ J l ^V.P)^ > (35) 
2 ^ j^L' 1 2 J 

The coefficients a' UL and are given by Eqs. ()82I83|) with the replacements k — > p and 
g — > p'. The same replacement is also required to calculate n 2 . 

Eq.(|34|l is the most general spin structure of the transition operator for a negative parity 
baryon B, consistent with symmetry principles. It also contains eight independent spin 
structures, but no central interaction is present in this case. 



n 
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Analogous to the positive parity case, Eq.fJMj) can be re-expressed in the form 



10 



A=l n,n'=0 



where the coefficients M^ n , are given by 



Mi 



Mi 



Mi 



M 



M 5 

nn 



Ml 



Ml. 



M 



M 9 

nn 



M. 
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Hi -\{u z + n 5 ) 

7^2 — r (^4 + ^4) 



^ln(l — <5n,o)^n',0 • 
p n (l - 5 n ,o)£n',0 , 

ni n 'S n>0 (l — 5 n f >0 ) 

Pn'^n,o(l — ^n',o) ) 



^2 t^ 3 ~ ~ ^n,o)(l ~ ^n',o) 



^7Pn^2n'(l ~ #n,o) (1 ~ ^n',o) , 



2 
1 

2 

-H8p' n h 2n >(l - 5 n ,o)(l - ^n'.o) 
-W8P^n 2n (l - 5 n ,o)(l - ^n',o) 



(36) 



(37) 



IV. APPLICATION: NEAR THRESHOLD AMPLITUDE IN pp M. 



As an application of the results of the previous section, we consider the reaction pp — > 
in the near-threshold energy region. We restrict to S and P waves in the final state, i.e., 
V = 0, 1. In contrast to the photoproduction reaction discussed in section II, here the Pauli 
principle restricts the initial pp state to (— ) 5+L = +1. 

For a positive parity baryon B, we then have six partial-wave states 



1 

On — ► On 



5 p , 1 p 3 p 3 p 3 77 3 p fQQ\ 



Here we use the notation 2S+1 L 



V x = A° , V 2 = (|)Alcos(0) 



2S ' +1 L'j. Then, the coefficients Pj in Eq.(|3U|) reduce to 



P 3 
V 7 



^) ™s(6) , V 8 = ) ^ , (39) 
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where the quantities Af,, Bf,, etc, are given in Appendix C. Here, again, the (p'/A) L ' 
dependence of the partial-wave matrix elements due to the centrifugal barrier is displayed 
explicitly. As in the photoproduction reaction discussed in section II, the scale A may be 
taken to be the four- momentum transfer at threshold, A 2 ~ — t = mErris' — rn 2 Nl where mg 
denotes the mass of the baryon B and thb', the mass of the positive parity baryon B'. 
For a negative parity baryon B, the possible partial- wave transitions are 

3 P - , 3 Pi - 3 Si , - 3 P , l D 2 - 3 P 2 . (40) 

With these partial- wave states, Eq.fSJ) reduces to 

Hi = , H 2 = A' Q n , H 3 = (?) B[ w sin(0) , H 4 = (?) C[ w cos(fl) , 



H 5 = , H 6 = C' 01 , H 7 = , H 8 = , (41) 

where the quantities A'f,' s , B'^,' s ', etc, are given in Appendix C. 

An interesting feature in Eq. (j39|) is that the amplitudes with spin-triplet initial states 
depend linearly on (p'/A) (i.e., the final state is in a P wave) if the parity of the baryon B 
is positive, whereas if the parity of B is negative (Eq.(|41jl). they do not depend on (p'/A) 
(i.e., the final state is an S state). This energy dependence is interchanged for amplitudes 
with spin-singlet initial states. The feature just mentioned is a direct consequence of the 
Pauli principle and total parity conservation, 



\S+L'+T 



Tl , (42) 



where the (— ) and (+) signs refer to the positive and negative parity baryon B, respectively. 
T denotes the total isospin (T = 1 for pp). In fact, as pointed out in Ref.|l5|. it follows 
immediately from Eq. (j42j) that the energy dependence of a partial- wave amplitude with 
spin-triplet initial state (S = 1) is given by an odd power of (p'/A) if the parity of B is 
positive since in this case the final state must be in an odd orbital angular momentum L', 
whereas it is given by an even power of (p'/A) if the parity is negative since V is even in this 
case. For completeness, we note that for the pn — > B'B reaction, this energy dependence is 
interchanged [lq . 

With the coefficients T>i and Hi given by Eqs. (|H9|41|) . it is straightforward to calculate any 
observable of interest in the pp — > B'B reaction using Eqs. f)3 1I32I36I37|) and the method given 
in Appendix D. Alternatively, one can calculate the observables using the matrix elements 
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given by Eq. (j24j) . in terms of which, some of the observables are also given in Appendix D. 
In particular, the cross section with the spin-triplet initial state Vs, as defined by Eq. (J§%|) 
in Appendix D, can be expressed as 

for a positive parity baryon B and, 

for a negative parity baryon B. Note that we have divided out the cross section by a factor 
of (p'/A) due to the final state phase space which introduces an extra p' dependence into 
the cross section. In the above equations 

Po = 2\B\\ 2 + 1 -{\C\\ 2 + \E{\ 2 ) -\U\\ 2 , 
A = -Po + l\A{\ 2 + l\Dl + El\*-2\Wl\ 2 , 

& = 2|A' 11 | 2 + |Cn 2 , (45) 

where U\ and Wl are defined in Appendix C. 

Eqs. (|43|44j) show that, apart from the p' dependence due to phase space, the spin-triplet 
cross section scales quadratically in p' (or equivalently, linearly in excess energy Q since 
y/Q oc p'), if the parity of B is positive, whereas it is constant if the parity of B is negative. 
This is the principal result of Ref . [lsj] . where the reaction NN — > YQ + has been investigated. 
Furthermore, the spin-triplet cross section exhibits a cos 2 (8) angular dependence in the case 
of a positive parity B, whereas it is isotropic in the case of a negative parity B. Therefore, 
an observation of a strong cos 2 (8) dependence in the measured angular distribution near 
threshold would imply a positive parity baryon B. An isotropic angular distribution, on the 
other hand, would be inconclusive about the parity of B. 

The spin-triplet cross section can be related to the spin correlation coefficients An as (see 
Appendix D) 

d( 3 ay,) 1 da , , . „ 

dSl = 45sj' 2 + + ' (46) 

with da/dVt denoting the unpolarized cross section and, An, the spin correlation coefficient 
with the spin orientations along the z-axis. Therefore, in general, Vs can be extracted 
experimentally by measuring the spin correlation coefficients A xx and A yy in conjunction 
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with the unpolarized cross section. Furthermore, it is immediate from Eg. (1951) that, at 

n 

threshold, as pointed out by Rekalo and Tomasi-Gustafsson in Ref. [16( in connection to the 
reaction pp — > E + + , A xx = A yy = — 1 for a positive parity baryon B while A xx = A yy > 
for a negative parity B. 



V. THE REACTION N + N^M + B' + N 

We now focus on the process N + N — > M + B' + N, where M stands for a spinless 
(scalar or pseudoscalar) meson and, B' for a spin-1/2 and positive-parity baryon. We start 
by making a partial-wave expansion of the corresponding reaction amplitude 

< S'M s ,\M(q,p';p)\SMs >= £ i L - L '~ l {l mi J' Mj,\JMj)(S' M S ,L' M L ,\J' M r ) 

x {SM s LM L \JMj)M^ SJ {q,p'-p) 
x Y lmi (q)Y L , MLl (p')Y* ML (p) , (47) 

where S, L, J stand for the total spin, total orbital angular momentum and the total angular 
momentum, respectively, of the initial NN state. Ms, Ml and Mj denote the corresponding 
projection quantum numbers. The primed quantities stand for the corresponding quantum 
numbers of the final B'N state. / and m/ denote the orbital angular momentum of the 
emitted meson and its projection, respectively, relative to the center-of-mass of the final 
two-baryon system. The summation runs over all quantum numbers not specified in the 
l.h.s. of Eq. (}4Tj) . p and p ' denote the relative momenta of the two baryons in the initial and 
final states, respectively, q denotes the momentum of the emitted meson with respect to the 
center-of-mass of the two baryons in the final state. We note that, in Eq.()47|). apart from 
the restrictions on the quantum numbers encoded in the geometrical factors, total parity 
conservation imposes that (— ) l+L+L ' = — 1 in the case of pseudoscalar meson production 
and (— y+ L + L ' = +1 in the case of scalar meson production. 

Eq. (j47|) can be inverted to solve for the partial-wave matrix element M^ / ,^' SJ (g,p / ; p). We 
have 

M ; i7 L ' SJ (g, p'; p) = t L ' +l - L (lmiJ'M. r | JMj) (S'M S ,L'M L , | J'M. r ) (SM s L0\ JMj) 

x J d(caB(e q ))Y? m (e q ,0) < S'M S r\M(q,p';i?)\SMs > , (48) 
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where, without loss of generality, the z-axis is chosen along p and q in the xz-plane; cos(0 9 ) = 
q-p. The summation is over all quantum numbers not specified in the l.h.s. of the equation. 
The most general spin structure of the transition operator can be extracted from Eq. (J47j) 

as 

M(q,p';p)= \S'M S , >< S'M s ,\M(q,p';p)\SM s >< SM S \ . (49) 

S'SM S M S , 

Inserting Eq. (}4Tj) into Eq. (}4T?]) and re-coupling gives 

^ \l' J' s'\ is' p j) 

x { }{ } ® B s]° ■ \ X HL')P lT , (50) 

<*P [J P 1 J [L S a) 

where [X(n,>)p z,] a is defined as 

[x mp it = [\ym ® Y L ,(p')r ® Y L (p)] a , (si) 

and contains all the information on the angular dependence of the transition operator. 

We now expand [Agi <8> Bs] a , for each tensor rank a, in terms of the complete set of 
available spin operators in the problem. Since the meson M produced in the final state is a 
spinless meson, the expansion of [As' <8> Bs] a is exactly the same as that for the N + N —* 
B' + B reaction discussed in section III and is given by Eq.([28p. 

Inserting Eq.(J2EJ) into Eq.([5U|) yields 

M(q,p';p) = H P S=0 + K X P S=1 + TZ 2 ■ {a x + a 2 ) + TZ 3 ■ {a x - a 2 )P s=0 + 7? 4 ■ (^ - a 2 )P s=1 
+ K\ ■ <g> a 2 f . (52) 

The above result is the most general spin structure of the spinless meson production operator 
in iViV collisions consistent with symmetry principles. The first two terms are the central 
spin singlet and triplet interactions, respectively. The third term is a tensor of rank 1 
and corresponds to the usual spin-orbit interaction. The fourth and fifth terms are also 
tensors of rank 1 but they describe the spin singlet— ^triplet and triplet— ►singlet transitions, 
respectively. The last term corresponds to the tensor interaction of rank 2. Apart from the 
fourth and fifth terms, all the other terms conserve the total spin in the transition. The 
coefficients multiplying the spin operators in the above equation contain the dynamics of 
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the reaction process and are given by 

n = E ii_iM (-) L W M ^ 0I, ^^^)[^')^] 



= ^E^ L '~'H J VW j^' J ' L l ^Mii:t\ q , v '- P )[x {lv)LL f 



. v , I' ./' 1 1 {3 J i r 



w)p l} 1 , 



^ = ~ E i L - L '- l H L - J+J ' + V'} [JfMtfFb, A p) 



X 

' 



v -, . L' J' 1 1 1 /? J , 

E^l , M [*<W l] 2 • (53) 



J (3 l\ LI 2 



Note that the summations over the quantum numbers in the above equation are restricted 
by total parity conservation and the Pauli principle: 



y+L+L' _ 



\L+S+T 



-1 , (54) 



where the ± signs refer to scalar (+) and pseudoscalar (— ) mesons, respectively; T denotes 
the total isospin of the two interacting nucleons in the initial state. If the baryon B in 
the final state is a nucleon, the summations in Eq.()53p are further restricted by the Pauli 
principle 

(-f+s'+r = _ x ; (55) 

where T' denotes the total isospin of the two interacting nucleons in the final state. 

The quantity [X(iL')p l}° m the above equation can be most easily evaluated by choosing 
the z-axis along the relative momentum p of the two nucleons in the initial state. We then 



21 



[X mp L ] a M = [ ^(PM a maM a )[YM) ® Y L ,(p')f 

a V47T 



have 



= [ ^}(PM a LO\aM a ) V {lm l L'M v \(5M a )Y lmi {q)Y U M Ll {p') . (56) 

The evaluation of [Xn L ,y L ] a can be further simplified if we choose the relative momentum 
p ' of the two nucleons in the final state in the xz-plane, in which case, Yjj m l , (p ') = 
H m l' y/(L'-M L ,y./(L> + M L ,)\[(2L' + l)/47r]V2p^ (p'.p). 

What we have done so far is completely general and applies to any spinless meson pro- 
duction in NN collisions. In the following, we treat the production of neutral and charged 
spinless mesons separately in the NN — ► MNN reaction. 



A. Neutral meson production amplitude in N + N^M + N + N. 

Since the only difference between the scalar and pseudoscalar meson production ampli- 
tude is that the sum I + L + V be even (scalar meson) or odd (pseudoscalar meson) as 
expressed in Eq. (I54jl . we restrict the following considerations to pseudoscalar meson pro- 
duction. The scalar meson production amplitude is easily obtained from the pseudoscalar 
meson production amplitude by changing the restriction on the orbital angular momenta 
from (-) l + L + L ' = -1 to (-y+ L + L ' = +1. 

In the case of neutral meson production, the available isospin operators in the problem 
are the usual isospin operators t% and r 2 together with the identity operator acting on the 
nucleon sector. We need to construct a scalar from these available operators. Since the total 
isospin of the two interacting nucleons is conserved in the production of a neutral meson, 
the isospin structure of the reaction amplitude is most simply expressed in terms of the 
isospin projection operator Pj- for the iso-singlet and iso-triplet transitions as T = and 1, 
respectively. Explicitly, Pt=o — (1 — T\ • r 2 )/4 and Pt=i — (3 + T\ • t 2 )/4:. 

We now note that, since the total isospin is conserved in the neutral meson produc- 
tion process, total parity conservation together with Pauli principle (Eqs. (l54|55jl ) demands 
that (-) l + s + s ' = -1 in the case of pseudoscalar meson production. This implies that the 
coefficients TZ^ and 7I4 in Eq. (j53|) involve the summation over even I only while all other 
coefficients in Eq.()53|) involve summation over odd I only. In what follows we will use the 
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superscript (e) or (o) in the coefficients in Eq. (|53Jl to indicate the restricted summation to 
even or odd I, respectively. 

With the above considerations, the (pseudoscalar) neutral meson production operator for 
a transition with a given isospin T is given by 

M T (q,p';p) = {n^P s=0 + n { ^P s=1 + n { 2 o) ■ {a x + a 2 ) + Tlf ■ {a x - a 2 )P s=0 

+ TZ^ ■ (a, - a 2 )P s=1 + Kf ] ■ [a, ® a 2 f)p T . (57) 

The above result is the most general spin-isospin structure of the neutral pseudoscalar meson 
production operator in NN collisions consistent with symmetry principles. 
In the total isospin basis, the pp and pn states are expressed as 

\pp > = \T= 1,M T = 1 > , 

\pn > = -=[\T = 1,M T = > +\T = 0,M T = >] , (58) 
v2 

so that, for pp collisions, we have 

M ppM (q,p';p) = M T=1 (q,p';p) , (59) 

while for pn collisions, 

M pnM (q,p ';p) = \J2 M T (q,p';p) ■ (60) 

T=0,1 

The transition operator for neutral scalar meson production is also given by Eq.()57|h 
except that the superscripts (e) and (o) in the quantities appearing in Eq.([57J) should be 
interchanged. 



B. Charged meson production amplitude in A^ + A^^M + A^ + A^. 

For charged meson production, the available nucleon isospin operators are the same as 
those for neutral meson production as discussed in the previous subsection. In addition, 
there is the isospin creation operator ir^ in the meson sector [\lm >= tt^O >], where the 
subscript m stands for the charge of the produced meson (m = +1, 0, —1). Again, we form 
all possible scalars from these operators with if m appearing once in each term. The isospin 
structure of the transition operator is then of the form 7r • O, where O stands for an isospin 
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operator of rank 1 in the two nucleon isospin sector. Three independent structures are 
possible for O: 

T = T' = 1 : 6 = (fi + r 2 ) , 

T = 0, r = 1 : 6 = {n - f 2 ) P T=0 , 

T=1,T' = 0: 6 = (ri-r 2 )P T =i. (61) 

In the isospin singlet— ^triplet and triplet— ^singlet transitions, total parity conservation 
and the Pauli principle (Eqs. (|54|55jl ) demand that (—y+ s + s — i n the case of pseudoscalar 
meson production. This condition is just opposite to that for either the isospin singlet- singlet 
or triplet-triplet transitions as discussed in the previous subsection. Taking into account this 
restriction together with Eq.([6ip. we have for a given transition T — > T' 

M T >T(q,p';p) = {^ O) p s =o + n[ o) p s =i + n { 2 0) ■ {s x + ff 2 ) + ■ (a x - a 2 )p s=0 

+ Hf ■ {a x - a 2 )P s =i + Hf° ] ■ [a, ® a 2 } 2 }ir ■ (n + f 2 ) 

+ {n { e) p s=0 + n[ e) p s=1 + Hf ■ ft + a 2 ) + ni o) ■ {B x - a 2 )p s=0 

+ K { : ] ■ - a 2 )P s=1 + Kf ] ■ [a x ® a 2 ] 2 }vr ■ - r 2 )P T . (62) 

Note that the coefficients of the spin operators in the two-nucleon isospin non-flip term 
have different restrictions on the angular momentum of the emitted meson relative to the 
corresponding coefficients in the two-nucleon isospin flip term. The above result is the most 
general spin-isospin structure of the charged pseudoscalar meson production operator in 
NN collisions consistent with symmetry principles. The transition operator for a specific 
reaction can be trivially obtained from Eq. (J62j) using Eq. (J58|) . 

The transition operator for charged scalar meson production is also given by Eq.([62p. 
except that the superscripts (e) and (o) in the quantities appearing in Eq.(JH2J should be 
interchanged. 

VI. NEAR-THRESHOLD pp -» ppM REACTION. 

In any particle production reaction, the near-threshold energy region is of particular 
interest due to the limited number of relevant partial-wave amplitudes. In this section, 
as an example, we consider the pp — > ppM reaction near threshold, where the produced 
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pseudoscalar meson M is primarily in an s-wave state. Then, considering only I = 0, the 
transition operators of the previous sections take particularly simple forms. Indeed, Eq. (j59|) 
becomes 



M, 



ppM\Hi 



p';p) = ll^ ■ {a 1 - a 2 )P s=0 + Kf • - a 2 )P s=1 



(e) 



(63) 



where the isopin matrix element has been already taken. From Eq. (|5Hjl . the coefficients TZ^ 



and TZ^ reduce to 



1Z 



(e) 



rJ [L}M^ L L (q,p';p)[X ( , 



0L')L' l] 1 , 



(<0 



with 



0L')L' L\ 



ll 



1Z 



[LL'] 



^= J2i L ' L '[L'}M t:i L '(q,p';p)[X { oL')L' l} 1 



(64) 



(4vr 



L>(L> + 1 



-(L'lL0\ll)Pl,(p ■ p')h x + (L'0L0\l0)P L ,(p ■ p')p 



■ (65) 



In the above equation, hi = \{p x p ') x p\/\p x p '\. 

Restricting the final two proton states to S and P waves, Eq. (j63|) reduces further to, 



Mpp M (q, p';p) = i{ \fr p' - (3 3 (P' ~ 3p • p'p)} ■ {B\ - a 2 )P s =o - Af(cfi- ^2)^=1 \ , (66) 



where 



fa 
fa 



1 



2(4vr) 3 / 2 
1 

2(4vr) 3 / 2 
1 /:> 



M^(q,p';p) , 
M^{q : p';p) 



(67) 



2(4tt) 3 / 2 V 2" 

are proportional to the partial-wave matrix elements corresponding to 3 Pq — > 1 S'os, 1 S'o — * 
3 Po s ; an d X -D2 ~~ ► 3 -P2S transitions, respectively Eq.()66|) agrees with the structure 

derived in Ref. based directly on symmetry considerations. We also note that the higher 
partial-wave contributions in pp —>■ pprj reported in Ref. ^ have been calculated based on 
the present method. 

The transition operator given by Eqs. (j57l59|) can be re-expressed in the form 

6 3 

M ppM = Y,Y1 M n V n (l)cv(2) , (68) 

A=l n,n'=0 
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where 



M 1 , 

Tin 






M 2 nn ' 






M nn' 


= \ n 2 + +^4 )Jn(l - 


"n.OjCn'.O ) 


M 4 , 

nn' 




^n',o)^n,0 ) 


ML 


= E £ nn'k(T^3 ~ ^/ )fc(l " 


" <Wl)(l - #n',o) 


M 6 , 

nn' 


= ^ ) (l _ 5niQ)(1 _ 5<()) j 





(69) 



with 



mirrt2m, 

r mim 2 



x | — - — (<J n ,i - m x m 2 b n ^) + (1 - |mi|)(l - |m, 2 |)6 ni3 
= E (-) m ^ ) (l^il^ 2 |2-m) 



(mim 2 , . . . r s mi .... . s 

X |? (m 2 On.lV, 2 + mi<V2<V,l) - ^(1 - |m 2 |)(O n ,l + WMn,2)<V,3 

- ^=(1 - |mi|)5 ni3 (5 n / i i + im 2 5 n / i2 )| . (70) 

Any observable of interest can then be calculated either from Eqs.()68J) or from Eq. (j47|) 
as given in Appendix D. 

In near-threshold kinematics with only s-wave mesons, we see from Eqs.flHHEHI) that the 
non- vanishing coefficients are M ~?' ,5 . If, in addition, the final two nucleon states are 
restricted to S and P waves, we have 

ML' = ^l{(P2-P3)Pn+m i (p-p')-Pl]Pn}(l-5 n , )S nlfi , 

mL' = ~4{ (/?2 " ^ )p ' n ' + WP ' p,) ~ W«'K o(1 - ' 

ML' = -\Y, £ nn'h{W2-^)p' k +WS-p')+Pl]h}^ , (71) 

k 

where (3^8 are defined in Eq.(|6TJ). 

Inserting the above results into Eqs. ()90l91l92|) (or alternatively using Eq.()47p and 
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Eqs.dnSMia)) yields 

^ = |/3i| 2 + |/3 2 | 2 + (3cos 2 ^+ 1)|/3 3 | 2 + 2(3cos 2 ^ - l)3^[/5 2 /3*] , 

—A - 

do du 

dn xx = dn yy 

= |A| 2 - |&| 2 - (3cos 2 fl + 1)|/5 3 | 2 - 2(3 cos 2 - , ( 72 ) 

which are equivalent to the results derived directly from Eq.(jHEJ) in Ref . [lsj] . From the above 
results, the final state (NN) S-wave contribution can be isolated via the combination: 

iSd+^-^-W. (73) 

which is nothing other than a consequence of the Pauli principle and parity conservation as 
discussed in section IV. 



VII. SUMMARY. 



Based on the partial-wave expansion of the reaction amplitude, we have derived the most 
general spin structure of the transition operator for the reactions 7 + N — ► M + B and 
N + N — > B' + B. Also, we have derived the most general spin structure of the spinless 
meson production operator for the N + N—>M + B' + N reaction. The present method 
used to extract the spin structure of the transition operator is quite general and, in principle, 
can be applied to any reaction process in a systematic way. The advantage of this method 
is that it relates the coefficients multiplying each spin operator directly to the partial-wave 
matrix elements to any desired order of the corresponding expansion. 
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VIII. APPENDIX A 

In this appendix we will determine the coefficients Oz/l, b L > L and cl>l in Eq.© as well 
as a! L i L and b' VL in Eq. (fT3|) . 
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Taking the scalar product of the last equality in Eq.© with [q®q] 2 , [k®k] 2 and [k®q] 2 , 
respectively, we have 

3m = 2cll'l + (3 cos 2 9 — + 2 cos9cl'l 

3v = (3 cos 2 9 — l)a,L'L + ^vl + 2 cos9cl>l 

3w = 2 cos9a L > L + 2 cos9bi>L + -(3 + cos 2 9)cl>l , (74) 



where cos 9 = k ■ q and 



u = \Y L {k) ®Y L ,{q)} 2 -[q^qf 
v = \Y L (k) ®Y L ,(q)} 2 -[k®k\ 2 

w = [Y L {k)®Y LI {q)] 2 -[k®q] 2 . (75) 
In order to arrive at Eq.([74jh we have also made use of the results 

[q®qf- [q®q] 2 = jj 

[k ® k} 2 ■ [q <g> q] 2 = ^(3cos 2 fl-l) 
2 

[k®q] 2 ■ [q®qf = -cos 9 

[k ® q} 2 ■ [k ® q} 2 = 1(3 + cos 2 9) . (76) 
6 

Eq. (J74*|) can be readily inverted to yield 

a-L'L = — j— \2u + (1 + cos 2 9)v — 4 cos 9w] 
sin 9 

bi'L = — t- \( I + cos 2 9)u + 2v - 4 cos 9w] 
sin 9 

cl>l = \-2 cos 9(u + v) + (3 cos 2 9 + l)w] . (77) 

sin 

Choosing the quantization axis z along k, the quantities u, v, and w defined in Eq. (|75]) 
can be expressed without loss of generality as 



u = ^^[LL'](L0L'0\20)P L (k ■ q) 
v = ^^{LL'](L0L'0\20)P L ,(k ■ q) 



[LL>] 

w = 



V^^'OIOI/O) ^(L0L'M\2M)(101M\2M)(L'MlM\lmi)Y lmi (9,0) . (78) 
47r i W ^ 
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Similarly to what has been done above, taking the scalar product of the last equality in 
Eq. (fT3*j) with [k <8> n 2 ] 2 and [q ® n 2 ] 2 , respectively, we have 

2r = a' L , L + cos 0b' UL 

2t = cos 9a' L , L + b' L , L , (79) 

where 

r = [Y L (k)®Y L ,(q)} 2 ■ [k®h 2 ] 2 

t = \Y L (k)®Y L ,(q)} 2 ■ [q®n 2 f . (80) 
In order to arrive at Eq. (f79"|) . we have also made use of the results 

[k ® h 2 } 2 ■ [k ® n 2 ] 2 = - 

[k®h 2 ] 2 ■ [q®h 2 ] 2 = ^cos9 . (81) 
Eq. (fT9"j) can be readily inverted to yield 



2 / 

a L'L = —^(r-tcosl 
sin 6 



h 'vL = -i(*-rcos0) . (82) 
sin 6 



Choosing the quantization axis z along k, the quantities r and t defined in Eq. fJHUj) can 
be expressed as 

1 [LL>] 



—i 



4vr y/L'(L' + 1) 



(L0L'l\21)Pl,{k-q) 



= -i^J^^L'OlOl/O) ^ (-) M (L0L'M|2M)(L , Mlm|/m i ) 

x (lmll|2- M)^ mj (0,O) . (83) 



IX. APPENDIX B 

In this appendix we will give the explicit expression for the quantities Ay, By and Cy 
in Eo.(P|) and A' L ,, B' L , and C^, in Eq.flHl)- 
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Using Eqs. (|8|lip for those states specified in Eq.ffTHj). we find 

f 




M 22 - V2M 22 



2 : 
02 



M 22 +M 1 \ 2 



1 3 » 
4tt V 5 13 ' 



II 21 I 13 

M^ 5 - Ml 2 + 5 + 



l 13 



5 « 
2 MA 2 



5 23 

M^ 2 



Similarly, using Eqs. (jl5ll7j) for those states specified in Eq. (f27?]) . we find 



IK 



B' = — 




J(1 -4 A /-M 10 



_> 2 



- v/^M^ + 4 



'12 



X. APPENDIX C 



The quantities Af,, B^,, etc, in Eq.(|3^j) are given by 

f 



4° 



P 

x) * 



El 



An 



M ooo 



1 1 



r M no + 3 M m + 5 M ii2] ; 



4vr3 
f f 
47r4 
f 3 
4^2^ 



^i°i 1 , 



-M- 



112 
11 



1 /35 w112 
— W — M 
4tt V 22 



13 > 



f 

47T 



M 



112 
13 



-M lT ° 



M^ 1 - 



M- 



112 
11 
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For convenience we define 



wl = 



73 

47T 
47T 



112 
11 



112 
11 



Similarly, the quantities A'f,' s , B'f,' s , etc, in Eq. fl4T|) are given by 



A' 11 



4tt 2 V 2 01 ' 



P 



A ' 



°0 



1 1 

4^2 

1 1 

4^2 



M 



102 
12 



12 



— ^M 010 
47t2 Mo1 



(87) 



(88) 



XI. APPENDIX D 



In order to calculate the observables directly from the transition operators in 
Eqs. (l29J34J52j) . it is convenient to express them in the form 

3 



m=EE M n v„(iK,(2) 

A n,n'=0 

where <Jo(i) = 1, cr\{i) = cr x (i), etc., for i-th nucleon. 
The unpolarized cross section is then given by 

^ = -Tr[MM^\ 
dtt 4 L J 



(89) 



(90) 



For the analyzing power we have 



A, A' n'=0 



fc,n=l 



(91) 
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where e kn i denotes the Levi-Civita antisymmetric tensor. We note that, from parity conser- 
vation, A x = A z = for the two-body reaction NN — ► B'B; this result also holds for the 
three-body reaction NN — > MB'N in the coplanar geometry. 
The spin correlation coefficient Ay is given by 

3 

= J2{ m i M oo( M "jT + K(M X jT}+Yl 2 *l M L(M*ye knj + MUM*ye kni ] 

A, A' k,n=l 
3 

- E M^(M^,re nkt e n/k/j } . (92) 

k,k' ,n,n'=l 

Of course, any observable can also be expressed in terms of the spin-matrix elements 
given by Eqs. (j24l47|) . We have, e.g., 

% = \ E I < S'M' S \M\SM S > | 2 , (93) 

S,S',M S ,M' S 

'' a A x = ^ E S'M' s \M\S0>< S'M' S \M\1-1>* 



S,S',M' S 
\1+S 



+ (-y +b < S'M' S \M\11 >< S'M' s \M\S0 >* 



' h A v = E %l< S'M' s \M\S0>< S'M' S \M\1-1>* 



S,S',M' S 



+ (-) 1+s < S'M' S \M\11 >< S'M' s \M\S0 >*} , 
^A z = ~ E [l < S"M£|M|11 > | 2 - | < S"M£|M|1 - 1 > | 2 ] , (94) 



S',M' S 



^A xx = \ H 1+S ~ 2Ms < S'M' S \M\SM S >< S'M' S \M\S - M s >* , 

If 4 ™ = I E R 1+5 ~ Ms < S'M' S \M\SM S >< S'M' S \M\S - M s >* , 

S,S',M S! M's 

^A zz = I E (") m/S l < 5'M^|M|^M S > | 2 . (95) 

S,S',M S ,M' S 

As has been pointed out in Ref.^|, a given initial state spin contribution to the cross 
section can be isolated by measuring some spin observables. In particular, using the spin- 
projection operator P$ as given in section III in terms of the Pauli spin matrices, it is 
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immediate that 

^M. = ±Tr[(MP s = )(MPs= y) = -Tr[MP s= oM ] ] 
I da 

~ 4 ~dTL ~ ~ y ~ zz > ' 
= l T r[(MP s=1 )(MP s=i y} = ~Tr[MP 5=1 Mt] 
X d(j 

= + A xx + A yy + A zz ) , (96) 

where 2S+1 a denotes the (initial state) spin-singlet or -triplet cross section as S = or 1, 
respectively. 

Also, Eql931 can be inverted to solve for the spin cross sections, d( 2S+1 Om s ) / d{VL) = 
J2s',M' a I < S'M' S \M\SM S > | 2 /4. We obtain Q 

IT 

,1 A xx Ayy A zz ^ , 





d(V ) 


1 da 




dVL 


= Idn 




d( 3 a ) 


1 da 




dn 




d{ 3 a x ) 




1 do" 


dtt 


dn 


~ 2dtt 



(1 "I" ^IX "I" Ayy A zz ^j , 

(1 + ^U, (97) 

where we have made used of the relation a = 1 a Q + 3 a + 3 <7i + 3 cr_i. Note that Eq. (f9Tj) is 
consistent with Eq. fjDTj) as 1 a = 1 a and 3 a = 3 cr + 3 o"i + 3 o-_x. 

From Eq.(|9*7|). it is immediate that the spin-triplet cross section defined as 

d(V s ) _ d( 3 a ) 1 fd^) d( 3 a^) \ 



dn dn 2 v rffi 

is given by Eq.flHJ). Note that d( 3 ai)/dVt = d{ 3 a^{)/dVt for NN -> £'£ by symmetry. This 
holds also for iViV — > MB'N in the coplanar geometry or for the cross section integrated 
over the emission angle of the one of the three particles in the final state. 
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